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NONLINEAR  ADAPTIVE  FLIGHT  CONTROL  WITH  A  BACKSTEPPING  DESIGN 

APPROACH 

Marc  L.  Steinberg^  and  Anthony  B.  Page^ 

Bldg  2187,  CST  5,  Naval  Air  Systems  Command,  Patuxent  River,  MD  20670 


Abstract 

This  paper  examines  the  use  of  adaptive  backstepping 
for  multi-axis  control  of  a  high  performance  aircraft. 
The  control  law  is  demonstrated  on  a  6  Degree-of- 
Freedom  simulation  with  nonlinear  aerodynamic  and 
engine  models,  actuator  models  with  saturation,  and 
turbulence.  Simulation  results  are  demonstrated  for 
large  pitch-roll  maneuvers,  and  for  maneuvers  with 
failure  of  the  right  stabilator.  There  are  substantial 
differences  between  the  control  law  design  and 
simulation  models,  which  are  used  to  demonstrate 
some  robustness  aspects  of  this  control  law.  Actuator 
saturation  is  shown  to  be  a  considerable  problem  for 
this  type  of  controller.  However,  the  flexibility  of  the 
backstepping  design  provides  opportunities  for 
improvement.  In  particular,  the  Lyapunov  function  is 
modified  so  that  the  growth  of  integrated  error  and  the 
rate  of  change  of  parameter  growth  are  both  reduced 
when  the  surface  commands  are  growing  at  a  rate  that 
will  likely  saturate  the  actuators.  In  addition,  the 
deadzone  technique  from  robust  linear  adaptive  control 
is  applied  to  improve  robustness  to  turbulence. 

Introduction 

In  the  early  1990's,  adaptive  backstepping  was 
developed  by  Kanellakopoulos  et  af  as  a  way  of 
designing  stable  adaptive  control  laws  for  a  broad  range 
of  nonlinear  systems.  Adaptive  backstepping  is  an 
approach  that  combines  Lyapunov  stability  theory 
with  the  substantial  advances  made  in  nonlinear 
differential-geometric  control  theory"  in  the  recent  past. 
The  basic  concept  behind  backstepping  is  to  use  some 
states  as  virtual  controls  to  control  other  states. 
However,  this  initial  approach  required 
overparameterization,  and  yielded  high  order 
controllers  that  were  not  very  practical  for  most  flight 


control  problems.  More  recent  work  by  Krstc  et  aP*^ 
developed  a  tuning  function  approach  to  eliminate 
overparameterization.  The  tuning  function  approach 
adds  computational  complexity  to  the  controller,  but 
yields  a  controller  of  the  same  order  as  the  number  of 
unknown  parameters.  The  primary  benefit  of  this 
type  of  controller  would  seem  to  be  that  it  allows  a 
wide  array  of  nonlinearities  to  be  incorporated  in  the 
controller  design,  and  has  proven  nominal  stability  and 
convergence  of  error.  The  types  of  nonlinearities 
included  in  the  control  law  design  could  be  either 
nonlinearities  in  the  system  model  or  nonlinearities 
chosen  to  meet  the  complex  design  criteria  associated 
with  flight  control.  In  the  past,  such  nonlinearities 
have  not  been  very  successful  in  flight  control^,  but 
the  powerful  theoretical  tools  associated  with 
backstepping  may  make  such  designs  more  feasible. 
Another  potential  advantage  of  backstepping  is  that  it 
converses  very  quickly  (when  it  does  converge)  because 
it  does"  not  have  the  lags  associated  with  parameter 
identification  for  conventional  adaptive  control 
approaches.  As  a  result,  it  may  be  effective  in  dealing 
with  damage  and  failures.  However,  there  are  concerns 
about  the  robustness  of  backstepping  designs.  Even 
more  problematic  is  that  basic  backstepping  designs 
tend  to  generate  very  large  effector  commands.  This  is 
a  serious  problem  for  flight  control,  due  to  the 
importance  of  actuator  saturations  in  aircraft.  Luckily, 
the  flexibility  of  backstepping  design  seems  to  provide 
opportunities  to  mitigate  this  problem.  Freeman  & 
Kokotovic^  and  Zhao  &  Kanellakopoulos ,  for 
example,  provide  some  interesting  approaches  to  reduce 
the  magnitude  of  control  commands  generated  by  a 
backstepping  control  law. 

In  recent  years,  backstepping  has  been  applied  to  a 
number  of  research  problems  including  a  1  Degree-of- 
Freedom  wing  rock  problem^  electric  motor  control^' 


’  Engineer,  Senior  Member  AIAA 
^Engineer,  Member  AIAA 

This  paper  is  declared  a  work  of  the  U.S.  Government 


and  is  not  subject  to  copyright  protection  in  the  United  States 


19980810  085 

DTIC  QUALITY  Ei3?EGTlD  I 


CLEARED  FOR 
OPEN  POBltCATION 

JUN  3  1998 

PUBLIC  AFFAIRS  OFFICE 
NAVAL  AIR  SYSTEMS  COMMAND 


I 


engine  control^^,  robotics^^,  neural  network 
control^"^,  and  ship  control’^'*^.  Another  interesting 
approach  towards  nonlinear  direct  adaptive  flight 
control  is  that  of  Kim  &  Calise*^,  Leitner  et  al^^  and 
McFarland*^.  This  approach  uses  neural  networks,  but 
provides  theoretical  proofs  based  on  similar  tools. 
Singh  &  Steinberg^®  developed  a  backstepping  adaptive 
control  law  with  increased  robustness  through  the 
incorporation  of  integrated  error  in  the  Lyapunov 
function,  and  applied  it  to  a  simple  nonlinear  aircraft 
simulation.  The  aircraft  simulation  was  a  constant 
velocity  model  with  no  disturbances  or  actuator 
models.  The  main  contribution  of  this  paper  is  to 
demonstrate  this  adaptive  backstepping  approach  on  a 
much  more  complex  simulation  model,  and  modify  the 
previous  design  to  deal  with  turbulence  and  actuator 
saturations.  Dealing  with  actuator  saturations,  in 
particular,  provides  a  good  venue  to  demonstrate  the 
flexibility  of  backstepping  design.  The  Lyapunov 
function  from  ref.  20  is  modified  to  reduce  integrated 
error  growth  and  the  rate  of  change  of  parameter 
estimates  when  the  actuator  command  is  approaching 
saturation.  Also,  in  ref.  20, the  only  difference  between 
the  control  law  design  model  and  the  simulation  model 
was  that  the  lift  and  drag  effects  of  the  control  surfaces 
were  not  taken  into  account  in  the  design  model.  In 
this  paper,  there  are  many  more  differences  from  the 
simulation  model,  so  the  results  in  this  paper  examine 
the  robustness  of  this  control  law  to  a  wide  range  of 
errors  in  the  model.  Finally,  the  control  law  in  this 
paper  makes  use  of  the  deadzone  technique  from  linear 
robust  adaptive  controF^  to  improve  robustness  to 
turbulence. 


Aircraft  Simulation  Model 

The  aircraft  simulation  being  used  is  a  high 
performance  aircraft  with  2  engines,  2  stabilators,  2 
ailerons,  2  rudders,  2  leading  edge  flaps,  and  2  trailing 
edge  flaps.  The  simulation  uses  the  standard  equations 
of  motion  and  kinematic  relations”^ 


u  = 
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The  components  of  the  aerodynamic  forces 
(/\,Fy,F,  )  and  moments  (/,,,/n^,,n^;)are  calculated 

from  table  look-ups.  Gross  thrust,  T,  is  calculated 
from  the  following  equation: 

T  =  [!  +  «! or + 

where  a^,  ^2,  c,  and  k  are  constants,  Fj  is 
calculated  from  a  table  look-up,  and  is  lagged 
throttle  position.  The  throttle  model  is  a  first  order 
linear  system  with  a  variable  time  constant  and 
variable  rate  limit  based  on  the  value  of  Pjj..  The 
actuator  models  are  2nd  order  linear  systems  with  rate 
and  position  limits.  The  turbulence  model  is  the 
standard  Dry  den  Gust  model  from  MIL-STD-1797A^^. 


Control  Law  Design  Model 

For  purposes  of  design,  the  full  simulation  model 
would  yield  a  control  law  that  was  far  too  complex  to 
be  practically  implemented.  Also,  it  was  felt  to  be 
important  to  deliberately  have  some  major  differences 
between  the  design  and  simulation  model  in  order  to 
examine  the  robustness  of  the  control  law.  As  a 
result,  the  following  model  was  used 
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Some  of  the  key  simplifications  made  in  this  model 
are  constant  velocity  (a  separate  auto-throttle  will 
attempt  to  maintain  this),  no  lift  and  drag  effects  of  the 
control  surfaces,  and  none  of  the  higher  frequency 
dynamics  (e.g.,  actuators).  Also,  the  stabilators  and 
rudders  will  only  be  used  collectively,  and  the  effects  of 
flaps,  which  are  scheduled  with  Mach  and  angle-of- 
attack,  are  ignored. 


parameters,  is  a  subset  of  the  state  vector,  and  u  is 
the  vector  of  control  effector  commands  where 

y  =  (<p,a,py 
co  =  (p,q,rf 
ri  =  9 

=  ia,p,(t),e) 

^1  = 

We  will  next  define  the  error,  e  ,and  a  function  s  that 
combines  error  and  integrated  error 

s  =  e  +  KqX^ 

Xs=fc^ 


Backstepping  Control  Law  Design 
This  design  roughly  follows  that  of  ref.  20  and  for 
more  detail  the  reader  should  consult  that  reference. 
However,  there  are  several  differences  from  that 
controller.  First  of  all,  the  Lyapunov  function  is 
modified  so  that  integrated  error  growth  and  the  rate  of 
change  of  unknown  parameters  are  reduced  when  the 
actuators  seem  likely  to  saturate.  Secondly,  several 
scalar  constants  are  turned  into  constant  matrices. 
Third,  all  controller  parameters  are  frozen  if  the 
actuators  saturate,  and  the  actuator  error  continues  to 
increase.  Finally,  a  deadzone  is  used  to  improve  the 
robustness  of  the  system  to  disturbances.  The  latter  2 
changes  have  not  yet  been  justified  theoretically  for  the 
backstepping  design,  but  have  been  for  linear  adaptive 
control  designs[21]. 

Following  [20],  the  aircraft  equations  of  motion  can  be 
put  in  the  form 

j  (JC, )  +  )w,  +  (x,  )0} 

d)  =  V'o  W  +  V'l  )« 

7]  =  ^0  W  +  )« 


where  is  a  positive  definite  symetric  matrix  and 
is  the  output  of  a  command  generator  that  is  a 
linear,  stable  3rd  order  system.  The  function  is  a 

scalar  bounded,  twice  differentiable,  positive  function 
with  bounded  derivatives  that  normally  equals  1,  but 
will  be  decreased  to  reduce  the  growth  of  integrated 
error  when  the  actuators  are  approaching  saturation. 
Choice  of  will  be  discussed  in  greater  detail  below. 

Taking  the  derivative  of  S  yields 

i  =  )wj  +  5(  JCj  )((0^  +  di)  +  V 

(0-(0^+G) 

v  =  %-yc  +  Kof,e 

In  this  equation,  CO ^  represents  the  desired  value  of  the 

virtual  controls.  However,  since  these  are  states  and 
not  directly  controllable  effectors,  there  will  be  an  error 

d). 

We  will  next  choose  0)^  through  the  use  of  the 
following  Lyapunov  function 


In  this  form,  y  is  a.  vector  of  the  outputs  that  will  be 
controlled,  CO  is  the  virtual  control  vector  used  to 
control  y,  7]  is  the  vector  of  uncontrolled  states, 
Wj,W2,  and  are  the  vectors  of  unknown 


Uj  f,s  +  wlL^w^)l2 

VV,  =  W-  —  Wj 
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where  Lj  is  a  positive  definite  diagonal  matrix  and 
is  the  estimate  of  W- . 

Taking  the  derivative  of  t/j  yields 

U^  =  +5(x,)(dj  +  <yJ  +  v(xp  j,,/,)] 

+  wl^w^  +fys 

Q)j  can  then  be  chosen  as 

(i)j  =5''(x,)[-A:„5-0^(x,)w,  -vix^,y^,fJ] 

where  is  a  positive  definite  diagonal  matrix. 

As  a  result  of  that  choice  of  CO^ 

Ui  =  -(^1  i/c  -  /<•  +  L,w,  ] 

If  we  putti)^  in  the  form 

(Od=y^od-¥id^i  +Wid^2 
then, 

<»  =  V^0«  +  +  V'2>2  +  )W 

¥oa  =¥o-  ¥od  +  ¥u^r  +  V'2aW2 

¥la  =  ¥ld 
¥2a  =¥l-  ¥2d 

We  will  next  chose  a  2nd  Lyapunov  function  such  that 

c/2  =  t/,  +  +  wlL^w^  +  h’JLjW,,  +  xlK^^x^)!! 

where  ^33  ^  positive  definite 

diagonal  matrices  and  is  given  by 

=  /c® 

Taking  the  derivative  of 

U2  =  -iKnf,  -  l)s^s  +  +  LiW,  +  f.y/lcbj 

+  W2[L,^2  +  fc¥L^]  + 

+a)^f\B^s  +  \ff^„  +  Du  +  ATjjX^]  +  f^cb^o) 


We  will  next  chose  the  following  control  law 

u  =  b-\-B^s  -  (XI/^  -  Wod  +  V'l^i 

+  (V^l  -  ¥2d)'^2)-  ^22^  -  ^33^^) 

Wj  (0^5  + (VA,>) 

^2  =  fAWL^ 

where  is  chosen  so 

[D{x,w^)-D{x,wJ\  =  ^^{x,u)w^ 

As  a  result,  since  and  K22  are  postive  definite 

diagonal  matrices  and  is  always  positive,  then,  if 
is  chosen  correctly 

f>2 

Following  the  approach  of  ref  20,  it  can  be  shown 
that  the  for  the  nominal  system,  s  ,5),  and  e  all  tend 
to  zero  as  t 

is  made  up  of  2  components.  The  first  is  a  fuzzy 
logic  component  and  the  second  is  a  3rd  order  linear 
stable  system,  chosen  such  that  meets  the  above 

requirements  for  the  range  of  possible  inputs  from  the 
fuzzy  logic  system  .  Fuzzy  logic  was  used  because 

it  is  an  easy  to  design  way  of  building  a  function  to 
transition  from  1  to  the  minimum  value  of  .  The 

fuzzy  logic  rules  are  as  follows: 

1)  If  actuator  position  is  normal  and  actuator  rate  is 
normal,  then  f-^  is  normal  (i.e.,  1) 

2)  If  (actuator  position  is  large  or  actuator  rate  is  fast) 
and  s  is  increasing  rapidly,  then  is  medium 

3)  If  (actuactor  position  is  very  large  or  actuator  rate  is 
very  fast)  and  s  is  increasing,  then  f-^  is  medium 

4)  If  (actuactor  position  is  very  large  or  actuator  rate  is 
very  fast)  and  s  is  increasing  rapidly,  then  f  is  small 

5)  If  (actuator  position  is  near  saturation  or  actuator 
rate  is  near  saturation)  and  (s  is  increasing  or  increasing 
rapidly),  then  is  small 

In  addition,  if  the  actuator  saturates  and  filtered  actuator 
error  is  increasing  above  a  threshold  value,  then  is 

immediately  set  to  zero.  This  violates  the  above 
stability  proof,  but  is  necessary  to  prevent  departures 
in  some  situations. 
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Simulation  Results 

All  simulation  results  were  generated  by 
Matlab/Simulink  v.  4.2.1  using  the  RK45  integration 
routine  and  the  full  nonlinear  simulation  model  with 
actuator  and  engine  models. 

Fig.  1  demonstrates  the  performance  of  the  control  law 
making  a  combined  60  degree  roll  and  8  degree  angle- 
of-attack  change  while  attempting  to  maintain  zero 
sideslip  at  Mach  .5  and  45,000  ft.  altitude  with  no 
turbulence.  The  solid  lines  are  the  commanded  values 
and  the  dotted  lines  are  the  actual  values.  This  is  a 
particularly  difficult  flight  condition  for  the  controller 
because  the  low  dynamic  pressure  can  lead  to  large 
control  magnitudes  that  could  saturate  the  actuators. 
As  you  can  see,  the  roll  response  is  very  good.  The 
alpha  response  is  not  quite  as  good,  but  this  is  partly 
due  to  the  fact  that  the  auto-throttle  was  not  perfect  and 
the  velocity  was  fluctuating  during  the  maneuver. 
There  is  some  slight  saturation  of  the  mdder  actuator. 
However,  this  is  much  less  than  would  occur  with  the 
basic  control  law  of  ref.  20  without  the  actuator 
saturation  mitigation  approach  of  this  paper.  Fig.  2 
shows  the  desired  and  actual  values  of  the  virtual 
control  inputs  for  this  case.  Fig.  3  shows  what 
happens  for  this  same  scenario,  with  f^=l  at  all 

times.  In  this  case,  the  actuators  rapidly  saturate,  and 
the  aircraft  departs.  Of  course,  with  this  initial  control 
law,  smaller  gains  could  have  been  used.  However,  the 
use  of  smaller  gains  led  to  very  poor  tracking 
performance  throughout  the  envelope. 

Figure  4  shows  the  same  maneuver  at  the  same  flight 
condition  with  a  failure  of  the  right  stabilator  at  1.5 
seconds.  There  is  only  fairly  slight  degradation  of  the 
roll  and  sideslip  response  and  modest  degradation  of  the 
alpha  response.  The  reason  for  the  good  response  can 
be  seen  in  Fig.  5,  which  shows  the  change  in 
estimated  stab  effectiveness.  Due  to  the  lack  of  filters, 
the  parameter  converges  to  a  new  reduced  value  very 
rapidly  after  the  failure  at  1 .5  sec. 

Tables  1  through  3  provide  some  error  statistics  for  a 
1 80  degree  roll  with  an  8  degree  angle-of-attack  change 
followed  by  a  return  to  the  initial  roll  and  angle-of- 
attack  angles  while  attempting  to  maintain  zero 
sideslip  in  moderate  turbulence,  as  defined  in  MIL- 
STD-1797  [22].  The  controller  performs  acceptably  at 
all  flight  conditions.  As  mentioned  earlier,  the  low 
dynamic  pressure  condition  had  the  most  difficulties 
with  departures  due  to  actuator  saturation,  but  the  high 
dynamic  pressure  condition  has  the  worst  tracking 
errors. 


Despite  the  improvements  made  in  mitigating  the 
impact  of  actuator  saturationsthrough  the  function, 

the  control  law  still  cannot  be  used  for  very  fast 
maneuvers.  Fig.  6  shows  an  attempt  to  command  a 
much  faster  roll.  In  this  case,  the  actuators  saturate 
despite  the  use  of  and  cause  the  aircraft  to  depart 

Conclusions 

This  paper  demonstrates  an  adaptive  backstepping 
flight  control  law  on  a  complex  high  performance 
aircraft  simulation.  The  control  law  is  demonstrated  to 
have  very  fast  convergence  properties  and  to  provide 
good  performance  in  many  situations,  including  the 
loss  of  a  stabilator.  It  also  demonstrated  good 
robustness  properties  to  a  wide  range  of  modelling 
errors.  However,  more  work  needs  to  be  done  to  reduce 
the  magnitudes  of  actuator  commands,  while  still 
maintaining  good  performance.  Achieving  good 
performance  generally  required  limiting  the  commands 
to  the  controller  to  prevent  departure.  Some  of  the 
potential  ways  of  dealing  with  this  could  include 
better  control  allocation,  improvements  in  the  use  of 
the  function,  and  changes  in  the  Lyapunov  function 

to  reduce  the  effective  gain  of  the  system  or  to  modify 
the  gain  based  on  dynamic  pressure.  Luckily, 
backstepping  design  provides  considerable  flexibility  to 
alter  the  Lyapunov  function  in  such  ways.  However, 
it  is  usually  not  clear  what  effect  many  changes  will 
have  on  the  system  and  better  theoretical  tools  and 
further  simulation  will  be  required. 
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Average  Absolute  Error 

Maximum  Absolute  Error 

Phi  (deg.) 

.74 

1.9 

Alpha  (deg.) 

.67 

1.2 

Beta  (deg.) 

.43 

1.4 

Table  2  -  Mach  .9,  5K  ft  altitude 


alpha  and  alphac  (deg)  beta  (deg) 


Table  3  -  Mach  .5,  45K  altitude 


Average  Absolute  Error 

Maximum  Absolute  Error 

Phi  (deg.) 

0.71 

2.6 

Alpha  (deg.) 

0.64 

2.2 

Beta  (deg.) 

0.49 

2.4 

Figure  1  -  No  Failure  Case  with  Actuator  Saturation  Mitigation  Activated 
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Figure  3  -  No  Failure  Case  without  Actuator  Saturation  Mitigation 
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